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\ Abstract 



We present a first numerical evidence for the existence of a novel magnetic 
condensate proposed recently by one of the authors in SU(2) Yang-Mills theory. 
^ ' In our framework, the spontaneously generated color magnetic field identified 

with the Savvidy vacuum has the microscopic origin and is a consequence of the 
intrinsic dynamics of the Yang-Mills theory. It strongly suggests the Nielsen- 
Olesen instability of the Savvidy vacuum disappears and the stability is restored 
without the need of the Copenhagen vacuum. The implications to the Skyrme- 
Faddeev model are also discussed. These results are obtained through the first 
implementation of the Cho-Faddeev-Niemi decomposition of the Yang-Mills 
field on a lattice. 
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1 Introduction 



In the SU(2) Yang-Mills theory, Savvidy [1] has discovered according to the renormal- 
ization group equation that a non-perturbative vacuum with dynamically generated 
color magnetic field H has lower vacuum energy density than the perturbative vac- 
uum. This is possible only for the non-Abelian gauge theory with asymptotic freedom. 
Immediately after this discovery, however, Nielsen and Olesen [2] have pointed out 
that the effective potential V(H) of the color magnetic field H, when calculated ex- 
plicitly at one-loop level, develops a pure imaginary part; The real part of V(H) has 
an absolute minimum at H = Hq ^ away from H = and satisfies the renormaliza- 
tion group equation in agreement with the Savvidy argument, while the non- vanishing 
imaginary part also satisfies the renormalization group equation without the renor- 
malization scale \x dependence. The presence of the pure imaginary part implies 
that the Savvidy vacuum becomes unstable due to gluon-antigluon pair annihilation. 
Since the energy eigenvalue E n of the massless off- diagonal gluons with spin S = 1 
and S g — ±1, in the constant external magnetic field H z := H\2 is given by 



the Nielsen-Olesen instability is also understood as originating from the tachyon mode 
n = 0,S Z = —1, i.e., the lowest Landau level with antiparallel spin to the external 
magnetic field, 



In fact, Eq becomes pure imaginary in the low-energy region k\ < gH . 

On the other hand, it is well known that in QED without asymptotic freedom, 
the non-zero magnetic field does not lower the vacuum energy and hence no magnetic 
condensation is expected to occur. Incidentally, external electric field always desta- 
bilizes the vacuum by causing electron-positron pair creation in QED and gluon pair 
annihilation in Yang-Mills theory. Therefore, no spontaneous generation of electric 
field is expected in both Abelian and non-Abelian gauge theories. 

The Nielsen-Olesen instability of the Savvidy vacuum was derived based on the 
one-loop calculation of the effective potential. Therefore, some people consider it as 
indicating unreliability of the lowest-order loop calculation, i.e., artifact of the ap- 
proximation. However, there have been published a huge amount of papers dealing 
with the problem of the unstable modes since the Nielsen-Olesen paper, including the 
stabilization by higher order terms [3]. Moreover, the same problem exists also in 
the supersymmetric Yang-Mills theory in which the higher-order loop corrections are 
absent, because the covariantly constant background field strength is not supersym- 
metric [4]. 

A way to circumvent the instability of the Savvidy vacuum is to introduce the 
magnetic domains with a finite extension into the Yang-Mills vacuum, in each of 
which the tachyon mode does not appear as far as k 2 z > gH z . This resolution is called 
the Copenhagen vacuum. However, the Copenhagen vacuum breaks the Lorentz 
invariance and color invariance explicitly. This issue has been re-examined recently 
by Cho and his collaborators [5]. 

What type of vacuum is allowed and preferred in the Yang-Mills theory is an 
important question related to the physical picture of quark confinement. Can the 
instability be resolved even in the one-loop level by a new mechanism? 
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First, it is instructive to recall the assumptions taken in Nielsen and Olesen [2]. 

1. The color magnetic field H has a uniform magnitude \H\ in spacetime and a 
specific direction if 12 = H z (The direction is identified with the quantization 
axis of the off-diagonal gluon spin). 

2. A background gauge is taken as a gauge fixing condition. (Note that the back- 
ground gauge is exactly the same as the Maximal Abelian gauge [7] which has 
been adopted in recent investigations on quark confinement based on the dual 
superconductor picture [8].) 

3. The off-diagonal gluons are treated as massless throughout the analysis. 

Now we would like to remind you of the facts which have been obtained by the 
recent investigations on quark confinement since 1990: 

1. In the Maximal Abelian gauge, infrared Abelian dominance [9, 10] and mag- 
netic monopole dominance are observed, as first confirmed [11] in the numerical 
simulations on the lattice, see [12] for reviews. 

2. The off-diagonal gluons acquire the mass M which is much larger than the 
diagonal gluon mass [13,14]. The off-diagonal gluon mass M measured on a 
lattice is M = 1.2GeV. See Ref. [15-21] for analytical works. 

In the previous work [6], the stability of the Savvidy vacuum has been re-examined 
by taking into account these facts and a scenario of eliminating the Nielsen-Olesen 
instability has been proposed to recover the stability of the vacuum: A novel type 
of color magnetic condensation originating from magnetic monopoles can occur and 
provides the mass of off-diagonal gluons in the Yang-Mills theory. Moreover, a novel 
magnetic condensation removes the tachyon mode of the off-diagonal gluon and the 
Nielsen-Olesen instability of Savvidy vacuum disappears to restore the stability of 
the magnetic vacuum, if the magnetic condensation is sufficiently large. 

The dynamical mass generation for the off-diagonal gluons enables us to explain 
the infrared Abelian dominance and monopole dominance by way of a non-Abelian 
Stokes theorem. These are quite natural and consistent results for understanding 
quark confinement, since the condensation of magnetic monopoles is the key concept 
in the dual superconductor picture. Therefore, quark confinement can be compati- 
ble with the stability of the Savvidy vacuum without resorting to the Copenhagen 
vacuum. 

The above claims were confirmed at least to one-loop order in the continuum 
theory by calculating the effective potential [6]. As a technical device, we have applied 
the Cho-Faddeev-Niemi (CFN) [22,23] decomposition to SU(2) Yang-Mills theory to 
extract the magnetic monopole degrees of freedom explicitly from the non-Abelian 
gauge potential. 

The purpose of this paper is to go beyond the previous analytical calculations and 
to confirm some of the above claims by using the numerical simulations on a lattice. 
This paper is organized as follows. In section 2, we review the CFN decomposition 
which plays a crucial role in this paper and summarize the analytical results and some 
predictions obtained in the previous papers [6]. In section 3, we argue how the CFN 
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variables on a lattice are defined to perform the numerical simulations on a lattice. 
Our definition of the CFN decomposition on a lattice reproduces the expressions 
of the continuum formulation from the lattice counterparts, in the naive continuum 
limit of the lattice spacing going to zero. Moreover, we simulate the lattice Yang-Mills 
theory without breaking the global SU(2) symmetry respected by the CFN variable. 
This section constitutes a crucial step to discriminate our approach from the other 
works which are apparently similar to ours. In section 4, we present the first results 
of numerical simulations based on the lattice gauge theory using the lattice CFN 
variables set up in the previous section. The first numerical evidence is obtained for 
the existence of two types of vacuum condensates, which supports the recovery of 
the stability in the Savvidy vacuum as claimed in [6]. The final section is devoted to 
conclusion and discussion. 

In Appendix A, we show how the gauge invariance of the Yang-Mills theory is 
expressed in terms of the CFN variables. Then we discuss how the gauge fixing 
is performed to eliminate the gauge degrees of freedom, especially in the Maximal 
Abelian gauge. In Appendix B, we summarize the relationship between the gauge 
fixing on a lattice and the continuum limit. It is shown explicitly that the gauge fixing 
procedures on a lattice which is actually used in the numerical simulations reduce to 
those known in the continuum formulation [6,22,23] in the naive continuum limit. 

2 Results and predictions from analytical works 
2.1 CFN decomposition in the continuum 

We adopt the Cho-Faddeev-Niemi (CFN) decomposition for the non-Abelian gauge 
field [22-25]: By introducing a unit vector field n(x) with three components, i.e., 
n(x) ■ n(x) := n A (x)n A (x) = 1 (A = 1,2,3), the non-Abelian gauge field ^(x) in 
the SU(2) Yang-Mills theory is decomposed as 

v^o 

=5^0) = Cf,(x)n(x) + g~ l d^n(x) x n(x) +X M (z), (2.1) 

where we have used the notation: C M (x) := c /i (x)n(x), B M (x) := g _1 <9 M n(x) x n(x) 
and V^(x) := C^(x) + B At (x). By definition, C^(x) is parallel to n(x), while B M (x) is 
orthogonal to n(x). We require X M (x) to be orthogonal to n(x), i.e., n(x) ■ X M (x) = 
0. We call C M (x) the restricted potential, while X M (x) is called the gauge- covariant 
potential and B M (x) is called the non-Abelian magnetic potential. In the naive Abelian 
projection, C M (x) corresponds to the diagonal component, while X M (x) corresponds 
to the off-diagonal component, apart from the vanishing magnetic part B M (a;). 
Accordingly, the non-Abelian field strength J^,(x) is decomposed as 

:= d^K ~ d u ^ + gs/p x < = + + ^X, - + gX^ x X„ (2.2) 

where we have introduced the covariant derivative in the background field V M by 
_D M [V] = _D M :— dfj, + gV M x, and defined the two kinds of field strength: 

Ejb, =E IMU n, := <9 M c„ - d v c^ (2.3) 

IV =d^ v - ctB M + gB„ x !„. (2.4) 
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Due to the special definition of B^, the magnetic field strength is rewritten as 

= -gM„ xM u = -g-\d»n x d v n) = H^n, (2.5) 
Hya, := ~g l n ■ (d^n x d u n). (2.6) 

where we have used a fact that M.^ is parallel to n. Moreover, is shown to be 
locally closed and hence it can be exact locally. In other words, we can introduce the 
Abelian magnetic potential for H^: 

:= -g^n ■ (d^n x d u n) = d^h u - d u h^. (2.7) 

Thus we can introduce two kinds of Abelian potential c M and and the corresponding 
Abelian field strength = <9 M c„ — d v c^ and = d^h v — dyh^. We call c M the 
(Abelian) electric potential and the (Abelian) magnetic potential (partial duality), 
because represents the color magnetic field generated by magnetic monopoles [16]. 
The CFN decomposition is useful to extract the topological configurations explicitly, 
such as a magnetic monopole (of Wu-Yang type), one instanton (of BPST type), and 
multi-instantons (of Witten type). The gauge invariance of the Yang-Mills theory in 
terms of the CFN variable is discussed in Appendix A. 

2.2 Advantages of our method using the CFN decomposition 

We enumerate some advantages and and characteristics in our treatment of the mag- 
netic vacuum using the CFN decomposition. (Some of them have already been em- 
phasized by Cho [5].) 

1. In our approach using the CFN decomposition, the direction of the color mag- 
netic field M^x) = H lu/ (x)n(x) can be chosen arbitrary at every spacetime 
point x by using a unit vector n(x) indicating the color direction. The Lorentz 
symmetry and color (global gauge) symmetry are not broken by considering 
||H|| := ^Jm~W^ = g~ l ^(n ■ (<9 M n x d u n)) 2 . It is invariant also under the 
color reflection, n(x) — > —n(x). 

2. This formalism enables us to specify the physical origin of magnetic conden- 
sation as arising from the magnetic monopole through the relation, H^ix) := 
— g~ l n(x) ■ (dfj,n(x) x d„n(x)). This gives a microscopic description of the dy- 
namically generated color magnetic field H^ v which is not necessarily uniform 
in spacetime, in contrast to the Savvidy, Nielsen and Olesen. 

3. We can discuss the implications to the Skyrme-Faddeev model [26] which is 
supposed to be a low-energy effective theory of Yang-Mills theory. This model 
is expected to describe glueballs as knot solitons. 

4. The non-Abelian Wilson loop operator can be rewritten in terms of the CFN 
variables through the Diakonov-Petrov version of the non-Abelian Stokes theo- 
rem [17,27]. Hence we can separate the contribution from the magnetic variables 
in the Wilson loop average to examine the magnetic monopole dominance. 



4 



2.3 Predictions 



In the previous work [6] the following issues have been discussed. 

1. A novel type of vacuum condensation (B^ • B M ) > can occur in addition to 



the magnetic condensation (||H||) := (^(B^ x B„) 2 ^ > 0. Here (B p ■ B M > > 
0, is called the magnetic condensation of mass dimension two and (||H||) > 
represents the spontaneous or dynamical generation of color magnetic field 
corresponding to the Savvidy vacuum. They are caused by gluonic interactions 
due to magnetic monopole degrees of freedom which are extracted by the CFN 
decomposition and are expressed through n, i.e., 

g 2 M^ ■ B /t = (d p n)\ (2.8) 



||H|| := y/U^ ■ IHLjj, = y/(gM fl x B,) 2 = (0„n x ^n) 2 . (2.9) 

2. If a novel type of magnetic condensation occurs (B M • B M ) > 0, then the off- 
diagonal gluons X p acquire their mass Mx through the relationship M\ = 
g 2 (B M • B M ) > 0. Then the infrared Abelian dominance and the magnetic 
monopole dominance follows immediately from this fact, supporting the dual 
superconductor picture for quark confinement. 

3. The energy level (spectrum) of the off-diagonal gluons is shifted by M x , i.e., 



E — > J k 2 + M\ — qHq. If the off-diagonal gluon mass Mx obtained in this 
way is sufficiently large so that 

M\ =/(B,-B,)>( 5 ||e||), (2.10) 

the tachyon mode is eliminated and the stability of the Savvidy vacuum is 
restored. Therefore, a criterion of stability restoration is given by 

M\ g 2 (B p -B„) (Q9 p n) 2 ) 
' g\\m~ g(\m\) ' U{d,n^d un f 



In fact, the above statements are supported from analytical works as follows. 
Even in the massive case, the existence of a magnetic condensation has been shown 
(gf||H||) > 0, based on the effective potential in the one-loop level (improved by the 
renormalization group) [6] where the Maximal Abelian gauge written in terms of the 
CFN variables, 

X := D^X, = 0, (2.12) 

is adopted. 

Then, the existence of another magnetic condensation, g 2 (B M • B M ) > 0, can be 
shown [6] based on a simple mathematical identity (B p x B„) • (B At x B^) = (B At ■ B p ) 2 — 
(B^ - B„) (B M - B„ ) , which yields a lower bound on B p -B At , (B M -B M ) 2 - (B M xM u )-(M^xM u ) = 
(B p -B v )(B /i -B l/ ) > 0, i.e., g 2 (M fl -M p ) > (g\\M\\), leading to a lower bound of the ratio 

r > 1. (2.13) 
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Then the tachyon mode is removed. But the possible zero mode can not be excluded 
by this bound. A stronger bound is obtained [6] by using the Faddeev-Niemi variable 
[28], g 2 (M /1 ■ Mfj) > y/2(g\\M\\) , which yields a better lower bound on the ratio, 

r > y/2. (2.14) 

This bound is also obtained by another method, see e.g. [29]. Thus the tachyon mode 
and the zero mode are removed. In fact, the effective potential V(||H||) is real- valued 
for r > 1. In particular, the r j limit reproduces the Nielsen-Olesen pure imaginary 
part, i.e., instability. 

For the above arguments to work in the rigorous sense, the existence of the mag- 
netic condensation ^(||H||) > must be shown in the full non-perturbative level 
beyond the loop calculation. This automatically leads to the existence of a novel 
magnetic condensation g 2 (B M • B M ), if such a vacuum is stable. The precise value of 
the ratio r is not yet determined. In fact, there is no theoretical upper bound on r, 
while the lower bound r > y/2 is known. Hence, we perform Monte Carlo simulations 
on a lattice to attack these issues. 



3 CFN decomposition on a lattice and n field en- 
semble 

We denote by the CFN- Yang-Mills theory the Yang-Mills theory written in terms 
of the CFN variables. The CFN- Yang-Mills theory has the local gauge symmetry 
SU{2)f ocal x [SU (2) / U {l)} e local larger than the original Yang-Mills theory, since we can 
rotate the CFN variable n(x) by angle (x) independently of the gauge transforma- 
tion parameter u)(x) of ^(x), see Appendix A and [30] for more details. In order to 
fix the whole local gauge symmetry, therefore, we must impose sufficient number of 
gauge fixing conditions. Recently, it has been clarified [30] how the CFN- Yang-Mills 
theory can be equivalent to the original Yang-Mills theory after the gauge fixing of 
the local gauge invariance in the continuum formulation. This idea is implemented 
on a lattice as follows. 

Now we discuss how to perform the CFN decomposition on a lattice and define 
the unit vector field n x to generate the ensemble of n-fields. In the whole of this 
paper, we restrict the gauge group to SU(2). 

3.1 LLG and new MAG 

First of all, we generate the configurations of SU(2) link variables {U x ^}, 

U X)fl = exp[-z'e#<c$^(x)], (3.1) 

using the standard Wilson action based on the heat bath method [31] where e is the 
lattice spacing and g is the coupling constant. 1 We use the continuum notation only 
for the Lie-algebra valued field variables, e.g., ^(x). 

1 It is possible to adopt different relationships between the link variable and the gauge potential, 
e.g., U x ,im — exp[— ieggf^x + je/t)]. However, the difference appears only in higher order terms in 
the lattice spacing e. They do not affect our main results and hence the difference is neglected in 
what follows. 
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Next, we introduce the functional, 



F LLG [U; Q] = 5>(1 - n U x>ti ) -\jd 4 x [«) Q (x)] 2 (e - 0), (3.2) 

where n U x ^ denotes the gauge-transformed link variable defined by n U X)fl := VL X U X ^VL X+ 
with a gauge group element Q x being an SU(2) matrix defined on a site x, and 
(■^fi) n (. x ) denotes the gauge-transformed potential defined by (s/y) Q (x) = Q(x)[£/ fJi (x)+ 
ig^d^ix) = ^(x) + D^[£/]u{x) + 0(u 2 ) for Q(x) = e iau w. Here the arrow in- 
dicates the naive continuum limit e — > of the lattice spacing e going to zero, see 
Appendix B. Then we impose the Lorentz-Landau gauge or Lattice Landau gauge 
(LLG) by minimizing the function Fllg[U] fi] with respect to the gauge transforma- 
tion Q x for the given link configurations {U Xj/M }, i.e., 

minimizingQi^G^; ^] minimizing^ J d A x [(^) n (x)] 2 . (3.3) 

In the continuum formulation, this is equivalent to imposing the gauge fixing condition 
d^^x) = 0. Thus this procedure determines a set of gauge-rotation matrices {fl x }. 
Note that the LLG fixes the local gauge symmetry SU(2)f ocal , while the LLG leaves 
the global symmetry SU(2)^ lobal intact. See Appendix A. 




Configuration update 



Figure 1: Lattice CFN decomposition obtained by imposing nMAG and LLG. 



Subsequently, we impose the new Maximal Abelian gauge 2 (nMAG) by minimizing 
the functional Fmag[U',G], U x ^ := n U x ^ defined by 

FmagP; G] = 5>(1 - a 3 G U X ^ a 3 G U^) -> / d'x [(^) G (*)] 2 , (3.4) 



2 This procedure is the same as the usual MAG. However, the meaning is totally different from 
the usual MAG, as shown in [301. 
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with respect to the gauge transformation G x , i.e., 

mmimizmg G F M ag[U',G] — > minimizing^ J d A x [(A^) G (x)] 2 , (3.5) 

where G(x) = e tge<yX \ Here the Cartan decomposition for srf^ix) has been used, 

U Xjfl = exp{-ieg^{x)} = exp{-ieg[a^x)T 3 + A«(x)T a ]}, (3.6) 

where a M is called the diagonal gauge field and A^a = 1,2) is called the off-diagonal 
gauge fields with the SU(2) generators T A = a A /2(A = 1, 2, 3). 

The nMAG breaks the local gauge symmetry [SU(2)/U(l)] e global and leaves the 
local U(l) 61 symmetry and the global U(l) e symmetry intact. The superscript 9 
indicates that this U(l) is not a subgroup of the SU(2) gauge group for the original 
Yang-Mills variable srf^. See Appendix A. Note that the nMAG breaks also the global 

symmetry [SU(2)/U(l)]%tai, while it does not break SU(2)^. 

The ensemble of n-fields is constructed as follows. See Fig. 1. It is shown in 
Appendix B that the minimization procedure of the nMAG leads to the construction 
of n according to 

n x := Gla 3 G x = n A a A , [n A = -tx[a A Gla 3 G x ]) . (3.7) 

This is because this MAG leads to the gauge fixing for the CFN variables as 

F MAG [U ] G]=Y, i < 1 - n ^ n -+M^ '=Fmag[UM - / d*x [(X M )(x)] 2 , (3.8) 

if we identify the link variable as 

U X)ll = exp{-ze^[C M (x) + E^x) + X M (s)]} = n U x ^, (3.9) 

which we call the lattice CFN decomposition. Here Fmag [U ] n] implies that MAG is 
also realized as the minimization with respect to n x . Even if the initial configurations 
U x ,n (or ^(x)) are the same, the CFN variables C M (x), B M (x), X^rr) are not neces- 
sarily the same if a different n x is adopted. Therefore, / d 4 x [(X M )(x)] 2 changes the 
value depending on the choice of n x or G x , although it has no explicit dependence on 
them. 3 

By imposing simultaneously the LLG and the nMAG in this way, we can com- 
pletely fix the whole local gauge invariance SU (2)f ocal x [SU (2)/[/(l)]f oca/ of the lattice 
CFN- Yang-Mills theory. The global symmetry SU(2) g j^ bal is unbroken. 

Here we distinguish two cases related to the global symmetry SU{2) g i oha i. 

3 A different interpretation is as follows. The different MAG functional is obtained for the CFN 
variable as 

F MAG [U-<^=Y J ^-n x U x ,»n x+ fil ll ) := F MAG [U;n] f d 4 x [(X M ) n (x)] 2 , (3.10) 

x,fj, J 

if we identify the link variable with the CFN decomposition. 

U Xtlx = exp{-ieg[C„(x) + M„(x) + X M (a;)]}. (3.11) 

Even in this case, the same gauge fixing condition is obtained, D^fVjX^ = 0, apart from the 
exceptional case u> = 0, see Appendix A. 
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3.1.1 SU (2) g i oba i-breakmg case 

If the numerical simulations are performed in such a way that LLG and MAG are 
close to each other [32], in the sense that the matrices G connecting LLG and MAG 
are on average close to the unit ones, i.e., — (A — 1,2,3), i.e., G x = G X I, for 
the parameterization of SU(2) matrices, 

G x = G° x I + iGia A , GlG^eR, £(G£) 2 = 1. (3.12) 

fi=0 

then we observe that n x = a 3 or = (0, 0, 1), namely, n x are aligned in the positive 
3-direction and hence the non-vanishing vacuum expectation value is observed as 

(n*) = M6 A3 . (3.13) 

This implies that the global SU(2) symmetry is broken explicitly to a global U(l), 
SU (2) global — * ^(l)giobai- In the two-point correlation functions, the exponential decay 
is observed for the parallel propagator 

(n 3 x nl) ~ (n 3 ) (n 3 ) + ce"™ 1 * 1 = M 2 + ce"™ 1 * 1 , (3.14) 

and for the perpendicular propagator 

^X>XWe- m 'N, (3 . 15) 

Z a=l 

but m and m! are slightly different, but nearly equal to 0.9GeV. This result was 
reported by [32] and confirmed also by our preliminary simulations [33]. 

3.1.2 S77(2) 9 ; ;, a rinvariant case 

Our main numerical simulations are performed as follows. In the continuum for- 
mulation, the CFN variables were introduced as a change of variables which does 
not break the global gauge symmetry SU{2) g i oha i or "color symmetry", which has a 
correspondence with the local gauge symmetry SU(2)i oca i in the original Yang-Mills 
theory. Hence the nMAG can be imposed in terms of the CFN variables without 
breaking the color symmetry. This is a crucial difference between the nMAG based 
on the CFN decomposition and the conventional MAG based on the ordinary Cartan 
decomposition which breaks the SU (2) g i b a i explicitly. See Appendix A. Therefore, 
we must perform the numerical simulations so as to preserve the color symmetry as 
much as possible. 4 This is in fact possible as follows. 

Remember that the MAG on a lattice is achieved by repeatedly performing the 
gauge transformations. In order to preserve the global SU(2) symmetry, we adopt 
a random gauge transformation only in the first sweep among the whole sweeps of 
gauge transformations in the standard iterative gauge fixing procedure for the MAG. 
This procedure moves an ensemble of unit vectors n x to a random ensemble of n x 

4 Whcthcr the color symmetry is spontaneously broken or not is another issue to be investigated 
separately. 
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which is far away from n x = (0,0,1), although this procedure might increase the 
functional Fmag- Then we search for the local minima around this configuration 
of n x by performing the successive gauge transformations. The first random gauge 
transformation as well as the subsequent gauge transformations are accumulated to 
obtain the gauge transformation matrix G by which n is constructed. Beginning with 
the LLG and ending with the MAG in this way, we can impose both LLG and MAG 
simultaneously. 

Our numerical simulations are performed by using the standard Wilson action and 
periodic boundary conditions under the following conditions. After the thermalization 
of 3000 sweeps starting with cold initial condition, we have obtained 50 samples of 
configurations at 100 sweep intervals. For LLG and MAG, we have used the over 
relaxation algorithm. 

The data of numerical simulations in Table 1 show the vanishing vacuum expec- 
tation value 

(n£) = (A =1,2,3). (3.16) 



Table 1: The magnetization < > on the 16 4 lattice at (3 — 2.4. 







Mean value 


Jack knife error(JKbin=2) 


< n 1 


> 


-0.0069695 


± 0.010294 


< n 2 


> 


0.011511 


± 0.015366 


< n 3 


> 


0.0014141 


± 0.013791 



Moreover, we have measured the two-point correlation functions defined by G A (x) 



A A 



no summation over A). The two-point correlation functions exhibit almost 
the same behavior in all the directions (A = 1, 2, 3), see Fig. 2. 
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Figure 2: The plots of two-point correlation functions 
lattice axis on the 16 lattice at [5 = 2.4. 



A A 



for A = 1,2,3 along the 



These results indicate that the global SU(2) symmetry (color symmetry) is unbro- 
ken in our main simulations, in contrast to [32]. This is a crucial point. 
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Table 2: Fitted parameters on the 16 4 lattice at /3 — 2.4. 





a(l) 


a(2) 


a(3) 


G 1 


0.0254 


0.2719 


-0.2462 


G 2 


0.0241 


0.2745 


-0.0213 


G 3 


0.0284 


0.2604 


-0.0271 



The data suggest the exponential decay, 

G A (x) = (n*i4) ~ ce" m|a;| (A = 1, 2, 3). (3.17) 

The exponential decay implies that there exists the mass gap in the theory. This is 
confirmed as follows. In Table 2, we have given the values of three fitting parameters 
a(l), a(2), a(3) when G 3 (x) is fitted to the cosh-function: G 3 (x) = a(l) * cosh(a(2) * 
(x — 8)) + a(3). Here a(2) corresponds to the mass gap. Since the physical scale 
e{(3 = 2.4) is 0.26784 in unit of the string tension a, the mass gap reads m p h ys = 

a{2)/e(j3) = 0.27/(0.26784 v ^ : ~ 1/2 ) = 1.008 x 440MeV = 440MeV. This should be 
compared with the SU(2) mass gap, M = 1.5GeV, which could be regarded as the 
lowest glueball mass [34]. 

3.2 Imposing LLG as preconditioning before MAG 

Finally, we explain why the LLG is imposed before taking the MAG. From the be- 
ginning, we could have imposed the MAG by minimizing the functional, 

Fmao[U; G'] = 5>(1 - a 3 G 'U^ a 3 G 'U^) - f <?x [(^) G '(x)] 2 , (3.18) 

with respect to the gauge transformation G' x , once the link variable configurations 
{U x ,fj,}, U XtlM = exp[— iegg/^x)} are generated using the Wilson action based on the 
heat bath method. This is equivalent to minimizing Fmag\U]ti\ with respect to n x : 

F MAG [U ] G'] = Y,^ 1 - n - U ^ n ^ U L) -FmagPM - f d'x [(X M )(x)] 2 , 

X,fJ, 

(3.19) 

where the following identifications are made: 

n x := G'ta 3 G' x = n*a A , (n* = Ui[a A G'la 3 G' x ]) , (3.20) 

and 

U Xtlt = exp{-ieg[C^x) + B M (z) + X^x)}}. (3.21) 

However, it is observed that the resulting ensemble of n x becomes random as 
characterized by the specific two-point correlation function 

= \ sABs **> Le -> = 0* + y) or l sAB ( x = v)> ( 3 - 22 ) 
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and the vanishing vacuum expectation value 

(ni) = 0. (3.23) 

There is no correlation among the field n x on the different sites. This is because the 
original link variables {U Xifl } are generated due to the gauge invariant original action 
and are distributed randomly along their gauge orbits. Therefore, the transformation 
matrix G' x becomes random in bringing the original gauge field configurations to the 
gauge fixing hypersurface. See Fig. 1. 

From the technical viewpoint, this difficulty is avoided if we begin with the ordered 
link variables {U X IM } by a preconditioning which eliminates the randomness. From this 
viewpoint, the LLG could be regarded as a preconditioning [32,36]. As we explained 
in the above, however, the LLG in our approach plays a more essential and a totally 
different role of specifying the CFN decomposition by combining LLG with the MAG, 
rather than merely removing the randomness, as emphasized in [30]. 

3.3 Discriminating our approach from the others 

Although the technique of constructing the unit vector field n x given above has 
already appeared, e.g., in [32,35,36], there is a crucial difference between our approach 
and others. In [32,35], the unit vector field n x was regarded as the field variable of 
the Skyrme-Faddeev model which is conjectured to be a low-energy effective theory 
of Yang-Mills theory. However, the precise relationship between the Skyrme-Faddeev 
model and the original Yang-Mills theory is still under debate. (The paper [32] 
concluded with the negative answer.) In contrast, our approach can identify the 
lattice field lattice version of the CFN field variable n(x) obtained by the 

CFN decomposition of the original gauge potential ^(x) in Yang-Mills theory. In 
fact, the naive continuum limit of the MAG on the lattice agrees with the MAG for 
the CFN variable [6], see Appendix B. To the best of our knowledge, such an explicit 
relationship has not been elucidated in the previous works including [32,35,36]. We do 
not assume any model written in terms of the unit vector field n x , which is regarded 
as an effective theory of Yang-Mills theory. 

In [32] , it is studied whether the identification of the Skyrme-Faddeev (or Faddeev- 
Niemi) model as a low-energy effective theory of Yang-Mills theory is efficient or not. 
The Skyrme-Faddeev model can have the same pattern of spontaneous symmetry 
breaking SU(2) — > U(l) as the nonlinear sigma model. Therefore, if such sponta- 
neous breaking of the global SU(2) symmetry occurs, two massless Nambu-Goldstone 
bosons appear and the mass gap disappears. This is because in the Skyrme-Faddeev 
model there are no gauge fields into which the massless Nambu-Goldstone bosons 
are absorbed through the Higgs mechanism. To avoid this unpleasant situation, the 
global SU(2) symmetry was explicitly broken in [32] by choosing the configuration in 
the neighborhood of n = (0, 0, 1) among a large number of local minima. This view- 
point is consistent with adopting the ensemble of n x aligned in a specific direction, 
since the n x in [32] is the field variable of describing the Skyrme-Faddeev model, 
which is not necessarily the CFN variable n x . On the contrary, the variable n x in 
our approach always denotes the CFN variable of the original Yang-Mills gauge field, 
without referring to the Skyrme-Faddeev model. This viewpoint does not lead to the 
immediate contradiction. The relationship of the Skyrme-Faddeev model and the 
Yang-Mills theory is discussed in the final section in our framework. 
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4 Numerical results: magnetic condensations 



We present the first numerical evidence for the existence of two vacuum condensates 
(#||H||) and (g 2 M^ ■ B M ), indicating the recovery of stability in the Savvidy vacuum. 

4.1 Setting up the simulations 

We define the lattice derivative [35] by 

A At n(x) := n(x + fi) — j(x)n(x) : = s^x), (4.1) 

which guarantees automatically the orthogonality condition n(x) ■ A M n(x) = n(x) ■ 
s(x) = on a lattice by choosing 7(2;) as 

j(x) — n(x) ■ n(x + pi). (4.2) 

This is not the case for the naive lattice derivative d^n(x) := n(x + fi) —n(x). Then 
B M (x) on a lattice is defined by 

gM^x) := s ft (x) x n(x) = A^n(x) x n(x) = n(x + /i) x n(x). (4.3) 

The squared g 2 W^(x) agrees with s^x) 2 just as in the continuum case: 

g 2 Ml(x) = (s,(x) x n(x)) 2 = 1 - ( 7 ,(x)) 2 = s,(x) 2 = (A,n(x)) 2 . (4.4) 

A simple calculation shows that 

n(x) x (s^(x) x s u (x)) = 0. (4.5) 

This implies that s^(x) x s u (x) is parallel to n(x) and does not have the components 
perpendicular to n(x). Therefore, it is natural to define gM^ u (x) and gH^ v (x) on a 
lattice by 

gM^x) := -(s M (x) x s v (x)) = gH^ v (x)n(x). (4.6) 

and 

gH^x) = gn(x) • H^(x) := -n(x) ■ (s^x) x s u (x)). (4.7) 

This implies the equality of the squared quantities: 

g 2 M 2 Jx) = (s» x s u (x)) 2 = [n(x) ■ (s» x s u (x))} 2 = g 2 H 2 ^(x). (4.8) 

Our numerical simulations are performed on the lattice with the lattice size 12 4 , 
24 4 , 36 4 by using the standard Wilson action for the gauge coupling (3 = 2.1 ~ 2.7 and 
periodic boundary conditions. Staring with cold initial condition and thermalizing 
50*100 sweeps, we have obtained 200 configurations (samples) for 12 4 , 36 4 lattice and 
500 samples for 24 4 lattice at intervals of 100 sweeps. For LLG and MAG, we have 
used the over relaxation algorithm. 
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0.5 



Figure 3: The Savvidy-like magnetic condensation: 
( x /(A M n x A„n) 2 \ versus j3 = 2.0 ~ 2.7 on 12 4 , 24 4 , 36 4 lattices. 



4.2 Savvidy-like magnetic condensation 

We have measured the magnetic condensation 

(</||H(a;)||) := (gy/U^x) ■ M^Jx)) = ^(A,n(x) x A v n(x)) 2 ) ■ 

by changing (3 on the lattices with different sizes. This corresponds to the Savvidy- 
like magnetic condensation, but it has microscopic origin written in terms of the field 
n(x) as a part of the gauge potential See Fig. 3 for the numerical values of the 
dimensionless magnetic condensation versus (3 = 2.0 ~ 2.7 on 12 4 , 24 4 , 36 4 lattices. 
We have also measured the squared magnetic condensation 

(g 2 \\n^\\ 2 ) ■= (9 2 ^v(x) -H^(x)) = ((A,n(x) x A,n(x)) 2 ). 

See Fig. 4. We can estimate the variance, (( p(x)|| - (||H(x)||)) 2 ) = (||H(x)|| 2 ) - 

(||H(a;)||) 2 , and the standard deviation a : = ^((||H(x)|| - (||M(x)||)) 2 ), as discussed 
in the effective potential. 

4.3 A novel magnetic condensation 

A novel magnetic condensation predicted in [6] 

(g 2 M,(x) ■ B M (x)) = ((A p n(x)) 2 ) . 

has been measured as shown in Fig. 5. The value is larger than the Savvidy-like 
magnetic condensation, as suggested by the analytical lower bound mentioned before. 

All data of two magnetic condensations are collected in Fig. 6 where they are mea- 
sured in units of the string tension by way of the lattice spacing e{(3) = y / cr( / 9)/(jphys, 
as a function (3 (Fig. 7) where <7 p hys = (440MeV) 2 is the physical string tension 
and &(/3) is the (dimensionless) lattice string tension (determined by the magnetic 
monopole part of the Abelian Wilson loop), see [37] for details. 
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Figure 4: The squared Savvidy-like magnetic condensation: (g 2 ||IHI|| 2 ) := (g 2 M^ lu • H M „) 
((A^n x A u n) 2 ) versus /3 = 2.0 ~ 2.7 on 12 4 ,24 4 ,36 4 lattices. 




Figure 5: A novel magnetic condensation ((/ 2 B At • B M ) := ((A p n) 2 ) versus /3 = 2.0 ~ 2.7 on 



12 4 ,24 4 ,36 4 lattices. 
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Figure 6: Two vacuum condensates in units of the physical string tension a p h ys versus the 
lattice spacing e on 12 4 ,24 4 ,36 4 lattices. 
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The numerical value M{j3) measured on the lattice for the quantity of mass di- 
mension one is translated into the physical value M p h ys through the relation M(f3) = 
M phys ■ e{/3), 



M, 



phys 



M{(5) M(/3) 



M{(3) 



The magnetic condensations of mass dimension two are translated as 

<(A p n(*)) 2 > 09) ((A p n(x)r) {(3) 



(4.9) 



(d p n(x))' 



phys 



""phys • 



(4.10) 



Both magnetic condensations of mass dimension two increase monotonically as 
the lattice spacing e decreases (or (3 increases). 



|3-£(lattice spacing ) 




Figure 7: Lattice spacing e versus j3 (inverse gauge coupling) in units of the physical string 
tension a p hys, reproduced from [37]. 



4.4 The ratio 

The precise ratio r := (g 2 M 2 ) / (g||IHI||) between two magnetic condensations is plotted 
in Fig. 8. Although the respective condensation changes considerably with decreasing 
in the lattice spacing e (or increasing j3), the ratio converges to a value r ~ 1.8 in the 
continuum limit e | 0. The obtained value of the ratio r ~ 1.8 supports the recovery 
of stability of the Savvidy vacuum according to the argument [6]. 

The fact that the increase of two vacuum condensates and the constancy of the 
ratio with respect to (3 suggests that the composite operators g 2 B 2 and <?||H|| besides 
the field n have non-zero anomalous dimensions which are nearly equal to each other. 
The anomalous dimension of the field n is obtained by calculating the correlation 
function (n A (x)n B (y)^ in the short distance \x — y\ or high energy-momentum region, 
just as obtained in the non-linear sigma model in two dimensions which has the 
asymptotic freedom [38,39]. The numerical determination of the anomalous dimension 
of the composite operator is possible in principle. However, this is still beyond the 
ability of our numerical calculations and to be reserved as a future problem. 
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Figure 8: The ratio r := (g 2 M 2 ) / (<7||H||) of two condensates versus the lattice spacing e on 
12 4 ,24 4 ,36 4 lattices. 



4.5 Lattice effective potential 

The probability distribution of the local operator is obtained by calculating the 
expectation value 



(%>-$(*))). 



(4.H) 



The effective potential is obtained from this distribution by taking the logarithm and 
changing the signature [40], 



(4.12) 



The effective constraint potential [41] is defined for the averaged operator $ 
V~ l J2xev $0*0 over the four-volume V by 



\ \ xev 



(4.13) 



The value of the composite field, at which the potential has a minimum or the field 
distribution is maximum, is equal to the value of the vacuum condensate. This 
argument can be easily extended to a number of operators, $i, $2; • • ■ and the effective 
potential V eff (<pi, <pz, - ■ ■ ). 

In our case, we can define two effective potentials written in terms of the values 
of two composite operators: 



V ef f (fo, to) = - In - g*B\x))6(</> H - g\\M(x 



(4.14) 



and 



V, 



eff y(p B 2,(p H 



) = - in (s - v~ l g 2 ®\x)) 5 L H -v-'Y.g 

\ \ xev / V xev 



x) 



(4.15) 
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Figure 9: The local composite operators 4>b^{x) '■= B 2 (x) and <Ph{x) ■= obey 
(left panel) the probability distribution independent of x: exp[— V e ff(4> B 2, 4>h)] , and (right 
panel) the effective potential V e ff(4> B 2,4>H), at (3 = 2.3 on 24 4 lattice for 500 samples. 








Figure 10: The averaged composite operators 4> B 2 := V~ l J2 x ev ~^ 2 ( x ) an d 4>H '■= 
V^ 1 J2xev l|H( x )ll obey (left panel) the probability distribution: exp[— V e ff{(j) B 2, < Aff)]) an d 
(right panel) the constraint effective potential V e ff((J) B 2,4>H), at /3 = 2.3 on 24 4 lattice for 
500 samples. 
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See Fig. 9 and Fig. 10 for the effective potentials obtained in the LLG and SU(2) 5 z 6 a /- 
invariant MAG. Our simulations have shown that the local potential (4.14) is indepen- 
dent of the point x and hence the spacetime average of the local potential is plotted 
in Fig. 9. 

The numerical calculations show that the support of Kf/(0B 2 ; </>#) and the dis- 
tribution are contained in the allowed region g 2 M 2 > g\\M.\\ and that the minimum of 
V e ff((f>B2,<j)H) an d the maximum of the distribution are indeed shifted from zero in 
the allowed region. These results clearly indicate the simultaneous existence of two 
vacuum condensates, although two operators B 2 and ||H|| are always greater than or 
equal to zero. 5 In the deconfinement phase, the minimum is expected to be at the 
zero value of the composite operators B 2 and ||H||. 

Thus the numerical results obtained in this paper confirm the qualitative result 
obtained by analytical calculations to the one-loop level in the previous paper [6]. 



4.6 Lorentz invariance on a lattice 



11=1 . 

11=3 i 
|i=4 



Figure 11: Check of Lorentz invariance. Plots of (left panel) V 1 J2xzv \-^x,n/ ^ or M = 
1,2,3,4 and A = 1,2,3, (right panel) V~ l J2xev (H x ,iu>) for \i,v = 1,2,3,4, on 24 4 lattice 
at (3 = 2.3 using 100, 300 and 500 samplings. 



The magnetic condensations measured so far are defined in the Lorentz invariant 
way from the beginning. On the lattice, the Lorentz invariance (Euclidean rotational 
invariance) is inevitably broken due to a non-zero lattice spacing. However, the 
rotation invariance by angle ir/2 exists even on the isotropic lattice. To see this 
discrete rotation invariance, we have measured the vacuum expectation values of a 
component of the Lorentz vector (b£^j and the Lorentz tensor (H Xjfll/ ) = (n x ■ M. Xiflu ) 
on a lattice, see Fig. 11. What they vanish is a necessary condition for the full 
Lorentz invariance of the continuum theory (b^{x)^ = and (H^ix)) = 0. The 

vacuum expectation values (^B^^j and (H X)IJ J) on a lattice should be zero reflecting the 
discrete rotation invariance. In fact, Fig. 11 indicates that the respective component is 

5 We can see that the distribution of ||H|| in Fig. 9 is consistent with the value of the stan- 
dard deviation a calculated from date of Fig. 4 according to er := ((||gH|| — (||<7H||)) 2 ) = 
V(||.gH|| 2 ) - (||.gH||)2, e.g., a = Vl-028 - 0.867 2 S 0.526 at p = 2.3. 
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Figure 12: The effect of lattice Gribov copies: r, (g 2 M 2 ) and (<7||H||) (from up to down) 
versus the number of how many times the random gauge transformations are performed for 
the original configurations, at (3=2.35 with the lattice size L = 8 4 for thermalization=3000, 
iteration=100, number of configurations=30. 



extremely small compared to the relevant vacuum condensates \B>^ M y and ^M XyfJiU 2 J 
denned in the Lorentz invariant (and global gauge invariant) way. Moreover, it is 
observed that the absolute value with the error are decreasing monotonically, as the 
number of samplings is increasing, as expected. 

Thus we conclude that the magnetic condensations and the dynamical generation 
of color magnetic field presented in this paper do not mean the violation of the Lorentz 
invariance in the Yang-Mills theory. This result is in sharp contrast with the original 
Savvidy and Copenhagen vacuum. 



4.7 Lattice Gribov copies 

In our calculations of magnetic condensations, we have also estimated the effect of 
lattice Gribov copies due to the program of performing the gauge fixing on a lattice 
[42,43]. We have used a standard iterative gauge fixing procedure for MAG and LLG. 
In such a case, gauge fixing sweeps may be stuck for some local minima of a gauge 
fixing functional. Different local minima give rise to different gauge transformations, 
but they can not be distinguished from the viewpoint of the iterative gauge fixing 
procedure. These are the lattice Gribov copies. To check the effect of copies to the 
magnetic condensations, we generate 30 of SU{2) configurations {U Xtfl } on 8 4 lattice 
at (3 = 2.35. Then, we generate 4 of gauge equivalent configurations (i.e., copies) via 
a random gauge transformation before performing the LLG. Using these gauge copies, 
we estimated the novel type of vacuum condensation ^B 2 ^, the magnetic condensation 

(||H||) and the index for the stability restoration r = ^ 2 B^ / (p||HI||). Fig. 12 shows 

the ratio r, ^g 2 B^,and (<7||H||). The horizontal axis show the number of times of a 
random gauge transformation. 

This result shows that the ratio is stable, although the respective condensation is 
a little affected by Gribov copies. Therefore, qualitative analyses given in this section 
will not be affected by Gribov copies. 
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4.8 Anatomy of dimension two condensates 



We define the spacetime average of the vacuum expectation value for the local oper- 
ator 0(x) by 



((O)):- 

Using the CFN variable, ^(x) - 



1 



d A x (0{x)) . 



V/,(x) 



X^x) 



squared potential is decomposed as 



[x 



[x 



(4.16) 



the 



; 2 » = ((^» + 2((V-X)) + ((Xj)) 
«B 2 »+2«B.X»- 



(K)) 



(«)■ 



(4.17) 



Following Zakharov et al. [44], the minimum of the squared potential ((^ 2 )) with 
respect to the gauge transformation is gauge invariant. It should be remarked that the 
operator X 2 is gauge invariant under the SU(2) local gauge transformation II, see [30]. 
Therefore, the difference — ((X 2 )) could have a gauge invariant meaning. 6 
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Figure 13: Anatomy of dimension two condensates. From top to down, positive ({e 2 g 2 £/ 2 }), 
((eVX 2 )), ((eVC 2 )), ((e 2 <? 2 B 2 )) , and negative ((eVX M -B M )). (Left panel) the bare values 
measured on a lattice, (Right panel) the ratio of the respective dimension two condensate 
{{& 2 )) to the original dimension two condensate ((s^ 2 }}- The presented values are raw 
data, i.e., bare values without any renormalizations. In the simulations, the gauge potential 
£^ii{x) is extracted from the link variable U x ^ using the linear definition, i.e, egsrf^x) = 
^(U x ,[i — U\ „). The simulations are performed on the 24 4 lattice under the same conditions 
as in the other simulations. 



Note that the cross term ((V ■ X)) = ((B ■ X)) is eliminated by a special choice of 
the gauge transformation II, u/ = u\\n (residual U(l) invariance) even after the nMA 
gauge is imposed, and hence («» - ((c 2 » + ((I 2 )) + ((X 2 )) = ((V 2 » + ((X 2 )). 
However, it does no longer hold in general after the complete gauge fixing. 

6 Recently, there appeared many papers [19,44-46,48-50] discussing the gauge invariance for the 
spacetime average of mass dimension two condensate (srf^) i.e., Among them, Slavnov 

[48-50] claims a stronger statement that ((^f )) can be gauge invariant and can have the same value 
independent of the gauge fixing adopted. 
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It is pointed out in [6] that the off-diagonal gluon condensation of mass dimension 
two ((X 2 )) 7^ leads to the Skyrme-Faddeev model as a low-energy effective theory 
of Yang-Mills theory. The existence was also assumed as a key ingredient in the 
studies [28,51]. 

In view of these, we have measured various dimension two condensates constructed 
from the CFN variables, including the off-diagonal gluon condensation ((X 2 )) ^ 0, 
for the first time, see Fig. 13. Here we have performed the CFN decomposition on 
a lattice according to (3.9). The numerical simulations show that the cross term 
((V ■ X)) = ((B • X)) does not vanish and becomes negative, ((B • X)) < 0, after the 
complete gauge fixing, i.e., the Landau gauge fixing in addition to the nMAG. In 
other words, X p is not in the mass eigenstate X M . However, the condensate ((X 2 )) = 
((X 2 )) 7^ has the same value. 

5 Conclusion and Discussion 

We have implemented the Cho-Faddeev-Niemi decomposition in the SU(2) Yang-Mills 
theory on a lattice. Performing the Monte Carlo simulation on a lattice based on this 
framework, we have obtained a first numerical evidence for the existence of a novel 
magnetic condensation ^B^\ in addition to another magnetic condensation (||H||) 
corresponding to the Savvidy-like magnetic field. We have confirmed the existence 
of the vacuum condensations by calculating the effective potential on a lattice and 
obtained the stable value for the ratio r ~ 1.8 in favor of stability restoration of the 
Savvidy vacuum according to the previous paper [6]. Moreover, it has been checked 
that the magnetic condensations in question do not break the Lorentz invariance. 

In the previous paper [6], we have argued that the stability of the Savvidy vacuum 
is restored due to the dynamical mass generation of off-diagonal gluons caused by a 
novel type of magnetic condensation (with mass dimension two) coming from mag- 
netic monopole degrees of freedom. The off-diagonal gluons acquire the dynamical 
mass, M\ = g 2 (B^y due to the existence of a novel magnetic condensation yS>fj 
and removes the tachyon mode of the off-diagonal gluon to cure the Nielsen-Olesen 
instability of the Savvidy vacuum, while the diagonal gluon remains massless. To 
really confirm this claim, we must check whether the off-diagonal gluon mass Mx 
determined by measuring the decay rate of the correlation function (^X^(x)X^ (y)"j 
agrees with the magnetic condensation g 2 (M 2 l ). In order to know the absolute value 
of the condensate ^B 2 ^}, we need to know more detailed behaviors of the propagator, 
e.g., the anomalous dimension of the field n. Analytical attempt of calculating the 
anomalous dimension is now in progress within the continuum formulation. 

The other vacuum condensation (X 2 ) is also important. In fact, the off-diagonal 
gluon condensation of mass dimension i? proposed in the MAG [19], (X p -X p ) = 2A 2 ^ 
0, in the present framework yields the mass term for the field B M or the kinetic term 
for n through the interaction term |B p • B At X p • X p : 

g 2 A 2 M ll ■ B M = A 2 (d ll n) 2 . (5.1) 

Therefore, the off-diagonal gluon condensation yields the Skyrme-Faddeev model [26], 
which has been proposed as a low-energy effective theory of Yang-Mills theory and is 
supposed to describe the glueball by the knot soliton solution. 
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A way to obtain a gauge invariant characterization of dual superconductivity in 
QCD is to calculate the Wilson loop average. This is in principle possible in the 
same framework using the CFN decomposition based on a version of the non-Abelian 
Stokes theorem [17,27]. 

It is also important to clarify the relationship between magnetic condensation 
discussed in this paper and magnetic monopole condensation as a source of dual 
superconductor, in order to confirm the magnetic monopole dominance. The issues 
are to be reported in subsequent papers. 

Remarks on the gauge invariance 

The numerical simulations performed in this paper are based on a new interpreta- 
tion of the CFN decomposition proposed in the paper [30]. As shown in [30], in order 
to determine the configurations of the n field, the complete gauge fixing is needed, 
and we have adopted the nMAG and the Landau gauge in this simulation. In this 
sense, the obtained configurations of n field depend on the gauge adopted. 

However, the following points should be remarked. 

1. We can choose an arbitrary gauge-fixing condition other than the Landau to fix 
the SU(2) local gauge symmetry II which remains after the nMAG, as clarified 
in [30]. 

2. The n field is not a directly measurable physical quantity, since it is a vector 
indicating the color direction at each spacetime point and is not a color singlet 
object. Therefore, even if the n field is subject to changes by taking different 
gauge fixing conditions, it does not cause the observable phenomena. Hence, 
this does not lead to any difficulty 

Therefore, the problem is whether or not the resulting changes of the n field 
influence the magnetic condensations in question. As already mentioned in ref. [6], 
we know that the operators B 2 and yHj^ whose expectation values are measured 
in this paper are not invariant under the full SU(2) local gauge transformation II, 
although they are invariant under the global gauge transformation II (color rotation). 
(After the nMAG, the theory has the SU(2) local gauge symmetry II, see [30]). In the 
operator level, the full SU(2) gauge invariant combinations are obtained by including 
the electric components as (E^ + H^) 2 or V 2 = (B M + C M ) 2 = B 2 + C 2 . Therefore, 
from the viewpoint of gauge invariance in the operator level, we should have measured 
the quantities (<J {E^ v + HL^) 2 ^ or (v 2 ^. Nevertheless, we have measured only (|B 2 ^) 

or {^H 2 ^ in this paper and avoided including the electric components. The reasons 
are as follows. 

1. The dimension two composite operators B 2 , ^H 2 ~ , and the dimension four 

operators EI 2 ^ are gauge invariant in the operator level under the local U(l)n 
gauge transformation, see Appendix A. 2. This is the same setting as the original 
approach of Nielsen-Olesen [2]. Moreover, they are also color singlets, i.e., 
invariant under the global SU(2) gauge transformation II (color rotation). 

2. In this paper we are interested in the magnetic contributions coming from the 
topological degrees of freedom expressed through the n field. 
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3. According to the conventional wisdom, the pure color electric field makes the 
vacuum unstable due to gluon-antigluon pair annihilations in gluodynamics, 
just as the pure electric field makes the QED vacuum unstable due to electron- 
positron pair creations. Therefore, the inclusion of the electric components 
and E M1 , under our identification could be other sources for the instability of the 
vacuum. This leads to difficulties in demonstrating our claim that the magnetic 
condensations stabilize the vacuum by eliminating the tachyon mode. 

4. Our simulations are inspired by the analytical calculation [6] to the one-loop 
order. In the pure magnetic case, the effective potential is obtained in the 
closed form. When the electric field and the magnetic field are simultaneously 
included, however, no one has succeeded to obtain the closed form for the ef- 
fective potential and the known expression is a cumbersome infinite series. 

5. After taking spacetime average and the vacuum expectation value or Yang-Mills 
average, the gauge variant operator in the strong sense could become gauge 
invariant in the week sense as discussed in [48-50]. This leaves a possibility of 
SU(2) gauge invariance for (^fy and {^Hlj^y 

6. Incidentally, the same type of a mathematical identity is applied to yield the 
lower bound on Vj. But the identity (V^ ■ V M ) 2 - (V^ x V„) • (V p x V„) = 
(V p ■ V„)(V /t • V„) > 0, yields not so beautiful and not so useful result: 

g 2 Wl = g 2 Ml + g 2 cl > g^ v + 2[c* {d u nf - c,c u d,n ■ d u n]. (5.2) 

Indeed, the effect of the electric field is important. We have succeeded to separate all 
the CFN variables n, c M and X M . Therefore, we can now calculate the relevant quanti- 
ties and estimate the desired electric contributions. Some of the results are presented 
in section 4.8 and at the workshop [52]. We plan to perform the detailed investiga- 
tion as the next work. However, we wish to avoid to present the details. Because, 
the inclusion of such materials makes the paper longer and the presentation could 
become rather incomplete. Without presenting the detailed numerical calculations, 
we can show that the other electric contributions increase the value of off-diagonal 
gluon mass and they are larger than the magnetic contributions, at least in the bare 
values before performing the renormalization by subtracting the perturbative ultravi- 
olet divergent part. Therefore, the inclusion of the electric part does not change the 
main claim of this paper: the tachyon mode is eliminated as a consequence of shifting 
the spectrum of the off-diagonal gluons upward due to the novel type of magnetic 
condensation (it is sufficiently achieved without the positive electric contribution). 

Acknowledgments 

The numerical simulations have been done on a supercomputer (NEC SX-5) at Re- 
search Center for Nuclear Physics (RCNP), Osaka University. This work is also 
supported in part by the Large Scale Simulation Program of High Energy Accelera- 
tor Research Organization (KEK). K.-I. K. is financially supported by Grant-in- Aid 
for Scientific Research (C) 14540243 from Japan Society for the Promotion of Sci- 
ence (JSPS), and in part by Grant-in- Aid for Scientific Research on Priority Areas 



24 



(B) 13135203 from the Ministry of Education, Culture, Sports, Science and Technol- 
ogy (MEXT). 

A Gauge invariance and fixing in the CFN variable 

A.l Gauge symmetry 

For the CFN decomposition, 

*fn(x) = Cp(x)n(x) + g^d^x) x n(x) + X A ,(s) (A.l) 
the restricted potential c M and gauge covariant potential X M are specified by to and 

c„(x) = n(x)-^(x), (A.2) 

X M (x) = g^n(x) x D^]n(x). (A.3) 
The second equation is obtained by making use of the fact that 

/yV]n := <9 M n + gV ^ x n = d^n + gB„ x n := D^Mjn = 0, (A.4) 
which yields 

Dfj,[s/]n = d^n + x n = d^n + gV M x to + gX M x to = gX M x to. (A. 5) 

Therefore, the gauge transformations <5c M , 5X M are uniquely determined, once the 
transformations Sn and Ssrf^ are specified. 

• The fact n(x) 2 = 1 urges us to consider the local rotation by an angle 0{x): 

Sn(x) = gn(x) x 0{x) = gn(x) x 6±(x), (A. 6) 

where 0±(x) are the perpendicular components of 0{x) with two independent 
components (n ■ = 0). For the parallel component 0\\(x) = 8\\(x)n(x), the 
vector field n(x) is invariant. Therefore, it is a redundant symmetry, which 
we call U(l)e symmetry, of the Yang-Mills theory written in terms of CFN 
variables, since c M (x) and X M (x) are also unchanged for a given stf^ix). This 
symmetry is the local SU(2)/U(1) symmetry and denoted by [SU{2) /U{l)] e local . 

• The invariance of the Lagrangian is guaranteed by the usual gauge transforma- 
tion: 

S^(x) = D^[sf\u{x). (A.7) 

This symmetry is the local SU(2) gauge symmetry and denoted by SU(2)f ocal . 

Note that u(x) and 6(x) are independent, since the original Yang-Mills Lagrangian 
is invariant irrespective of the choice of 0(x). 

For later convenience, we denote the above transformations by Sg and 8^: 

(1) ^ (to - = 0), ui = 0: 

5 e n(x) = gn(x) x 0(x), 8 e ^(x) = 0. (A.8) 

(2) = 0, u ^ 0: 

6 u n(x) = 0, 5^(x) = D^]uj(x). (A.9) 

Then the general gauge transformation of the CFN variables is obtained by combining 
So and S w . 
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A. 2 Local gauge transformations I and II 

In the papers [5,6], two local gauge transformations are introduced by decomposing 
the original gauge transformation, 5 u ^^{x) = D fl [s^]u}(x). 7 
Local gauge transformation I: 

S u n =0, (A.lOa) 

S u c„ =n ■ D„\sf\u», (A.lOb) 

(J w X M =D M [^]w-n(n-D M [^«), (A.lOc) 

=>6 u B lt = 0, 6 u V li = n{n-D li [s/\u). (A.lOd) 

Local gauge transformation II: 

=gn x a/, (A. 11a) 

fe=n-^'. ( A - llb ) 

<£X„ =gX^ x w', (A. 11c) 

=^B M = [B]u/ - (n ■ fyw>> = (A. lid) 

The gauge transformation for the field strength can be obtained in the similar way 
as follows. 

Local gauge transformation I: 

cLE^ =n6 u Ep, = n{d^n ■ D v [srf\u) - d u (n ■ D^]u>)}, (A.12) 

SJE.^ =fwJ w ^ = 0. (A. 13) 

Local gauge transformation II: 

=gE^ xu' + n{d„{n ■ d v J) - d u (n ■ d^')}, (A. 14) 

5'JB.p, =gU fll/ x u/ - n{d„(n ■ d„u>') - d v {n ■ d^')}. (A.15) 

This implies the transformation for the sum 

<^(E^ + H^) =g(E^ u + W^) x a/, (A.16) 

leading to the full SU(2)n invariance: 

^(E M , + H^) 2 =0. (A.17) 

Moreover, we can show that 

5' U (D^Y]X U - A,[V]X„) =g(D,[¥]X u - D U [Y]XJ x w', (A.18) 

c^(X M x X,) =g(X^ x X„) x u>', (A.19) 



7 The gauge transformation I was called the passive or quantum gauge transformation, while 
II was called the active or background gauge transformation. However, this classification is not 
necessarily independent, leading to sometimes confusing and misleading results. 
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which lead to the full SU(2)n gauge invariance: 

5' u (D li \y\X v - D v \V\Xrf = 0, 5'JX^ x X,) 2 = 0, (A.20) 
<^[(E^ + H M „) ■ (gX^ x X,)] = 0. (A.21) 

In particular, when u)'(x) is parallel to n, i.e., u/(x) = 6'(x)n(x), we obtain 
Local U(l) gauge transformation II for uj'(x) = 9'{x)n{x): 

5' g n =0, (A.22a) 

=d^9\ (A.22b) 

=#X M x 9'n, (A.22c) 

=^1,, = 0, 5' e Y^ = nd^'. (A.22d) 

Note that n and B M are invariant under the U(l)n gauge transformation II, while 
c M transforms as the U(l)n gauge field. It is easy to show the local U(l)n gauge 
invariance for the field strengths: 

= 0, S'fir = 0, (A.23) 
which is also consistent with the initial definitions: 

V = n(d M c - cU M ), U„ u = -gM^ x B„. (A.24) 



Therefore, the dimension two composite operators B 2 , yH 2 tJy , and the dimension 
four operators are gauge invariant under the local U(l)n gauge transformation. 



A. 3 MAG as a partial gauge fixing 

The gauge transformation I defined in the previous paper [6] is nothing but 8^. On the 
other hand, the gauge transformation II has been defined in [6] as a gauge transforma- 
tion such that it does not change X 2 . To see this, we consider the gauge transformation 
of X 2 . Since the relationship (A. 5) leads to 

X 2 = g~\n x D^]nf = g~ 2 {{D^]nf - (n • D^]n) 2 } = g~ 2 {D^]n)\ 

(A.25) 

the gauge transformation of X 2 is calculated as [30] 

5 l -Xl = g-\D^]n) ■ {n x D^O 1 - - «-*-)}, (A.26) 

where we have used (A. 8) and (A. 9). Therefore, it turns out that the gauge transfor- 
mation II corresponds to a special case (x) = a;- L (x). 
The average over the spacetime of (A.26) reads [30] 

5 J d A x l -Xl = J d A x{6 L - u 1 ) ■ D^]X^ (A.27) 
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where we have used (A. 3) and integration by parts. Hence the minimizing condition 

= 5 J d 4 x^X 2 t (A.28) 

for arbitrary d^- and uj 1 - yields the differential form: 

Fma = X ■= D^YjX^ = 0, (A.29) 

which reproduces exactly the MAG for the CFN variables [6]. Therefore, the min- 
imization condition (A.28) works as a gauge fixing condition except for the gauge 
transformation II, i.e., 0- L {x) = a;- L (x). 



B Lattice CFN variables and gauge fixing 

B.l Continuum 

We show that for the CFN decomposition, 

< = V M + X M = C M + B A4 + X /1 , C M = c M n, B tl = g-%nxn, (B.l) 
the equality holds, 

{D^]nf = g 2 Xl- (B.2) 
In other words, X 2 is rewritten in terms of srf and n. This is shown as follows. 

(D^]n) 2 = («9 M n + x nf 

= (<9 M n) 2 + 2gd il n ■ (^ x n) + g 2 (^ x n) • x n) 
= (<V*) 2 + 2g^ ■ (n x dp-n) + g 2 (^ ■ <)(n • n) - g 2 (^ ■ n) 2 
=g 2 Ml - 2g 2 {C ll + B M + X M ) • % + /(C M + % + X M ) 2 - /(Cg 2 
=g 2 Ml - 2g 2 Ml - 2#% ■ M„ + g 2 (M fl + X M ) 2 

=<? 2 X 2 , (B.3) 

where we have used (A x B) • (C x D) = (A ■ C)(B • D) - (A • D)(B • C). 

Another (simpler) way of showing the equivalence between (D^[£/]n) 2 and X 2 is 
as follows. By making use the fact that 

D^[¥]n := d^n + g Y^xn = d„n + gM^ x n := D„\B]n ee 0, (B.4) 

we find 

D^]n = d^n + gs/^ xn = d^n + gV^ xn + gX^ xn = gX^ x n. (B.5) 
This fact leads us to the equivalence, 

p M K]n) 2 = g 2 (X^ x n) 2 = g 2 (X, • X M )(n • n) - g 2 (X^ ■ n) 2 = g 2 Xl. (B.6) 
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Hence, X M is rewritten in terms of ^ and n 

X M = g- l n x D^]n. (B.7) 

This is also the case for c M , 

c M = n ■ rfp. (B.8) 

As shown in Appendix A, we can impose the gauge fixing condition by minimizing 
the following functional under the local gauge transformation: 

= 5^X1 (B.9) 

which leads to the differential form of the MAG condition in the CFN decomposition 

D M [V]X M = 0. (B.10) 
The MAG condition can also be derived by minimizing the functional 

o = 5/J(A*Kh) 2 - (B.ll) 

This is confirmed by explicit calculation and we leave it for the reader as an exercise. 
B.2 Lattice 

We show that the link variable on the lattice is identified with the CFN decomposition 
of the gauge potential as 

U Xtfl = exp{-ieg^(x)} = exp{-ze^[C M (a;) + M^x) + X M (x)]}. (B.12) 

In fact, we recover (D^[£/]n) 2 in the naive continuum limit from the lattice functional 

Fmag '■= — 

Y.u-nJ\,.,,n,.. i ,r-. (B.13) 
In fact, by expanding the exponential U Xtfl = e- ie ^( x ) into the Taylor series, we obtain 

=y - \e 2 {{d^nf - 2g{d^n xn)-< + g 2 (n ■ n)(^ ■ - g 2 {n • ^) 2 } + 0(e 3 ), 



2 4 



eV^ + 0(6 3 ), (B.14) 



where the summation over /j, = 1, • • • , D should be understood and the order e terms 
cancel. Therefore, we can obtain the MAG in the CFN decomposition by minimizing 
the functional Fmag with respect to the gauge transformation under the identification 
(B.12). 
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In particular, the naive MAG for the usual Cartan decomposition is obtained from 
minimizing the functional 



D 



-i e y{«-o-(^-o 2 } + o(6 3 ) 
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=§ -\c 2 9 2 a;a; + o(6% (B.i5) 

for the link variable 

9 U Xilx = ex V {Tieg^(x)} = exp{Tie<?[a M (x)T 3 + A*{x)T a ]}. (B.16) 

In the above calculation of the naive continuum limit, we have used the following 
formulae for the trace of the product of the generators in the SU(N) algebra. 

tv(T A T B ) =\s AB , (B.17) 

tr(T A T B T c ) =- A (if ABC + d ABC ), (B.18) 

tr(T A T B T c T D ) =—5 AB 5 CD + \if ABE + d ABE )(if CDE + d CDE ). (B.19) 
AN 8 

They are obtained by the repeated use of 

T A T B = ±[T A ,T B } + 1 -{T A ,T B } = l -f*BD T D + 1 n 6 AB + jABD^ 2Q) 

For SU(2), they are simplified as 

tr(T A T B ) =-6 AB , (B.21) 
ti(T A T B T c ) =\ie ABC , (B.22) 

tr(T A T B T c T D ) =-5 AB 5 CD — \ e ABE e CDE — Ia^aC- _ Ia^Ca- 8 - _i_ Ia ad a bc< 
1 ' ~8 8 ~ 8 "8 8 

(B.23) 

They are obtained by the repeated use of 

T A T B = * ABD T D + \ 5 AB_ (R24) 

2 4 K J 

In order to obtain this result (B.14), we must symmetrize the expression, i.e., 

Fmag ■= - | ^tr[nJJ x ^n x+li Ul^ + n x - fJ ,U x -^ IJ ,'n x Ul_ llifl ] 

- tr[n x U x ^n x+ ^U ] X41 + n x U x - fl ,n x - l) ,Ul_ ll \, (B.25) 
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since +/i and —pi should be treated on the equal footing. Then the kinetic term for 
n is derived as 

]- tr [n x n x+fl + n x ^^n x ] =- $^tr[n x (n a .+ |1 + n,^)] 

x,i± x,fj, 

= Wt4n x (2n x + e 2 d^n x )} 
1 1 



E 



X 

£> 1 



(B.26) 
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where we have defined the forward derivative and the backward derivative by 

&${x) := [<f>{x + »)- 0(x)]/e, &^{x) := [<f>(x) - <f>{x - /i)]/e, (B.27) 
and the integration by parts, 

£ f(x)&lg(x) = -£ fl^/(x)y(x). (B.28) 

The lattice d'Alembertian is defined by 

A(x, y) := -<9^5(x -y) = -&^5(x - y). (B.29) 
The other terms can be calculated in the similar way. 

B.3 Remarks 

We show that both (D^[£/]n) 2 and are invariant 

5 II (D^}n) 2 = g 2 8 n Xl = 0, (B.30) 

under the local gauge transformation II of the CFN variables: 

n(x) — > U(x)n(x)W(x) := n'(x), 
X^x) - U{x)X ti (x)U^(x) := X'^x), 
D^\[x) -> U{x)D^]{x)U\x) := D^]'(x) } (B.31) 

where U(x) = exp(igu'(x)). Note that n(x) and X^(x) transform in the adjoint 
transformation under the gauge transformation II. In particular, D^[£/](x) transforms 
in the same way under the gauge transformation I and II, since it is written in terms 
of the original variable ^ which transforms as 

sf^x) -> U{x)^U\x) + -U{x)8^U\x) := ^(x). (B.32) 
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The signature in front of stf is important, since 

(D^-^n) 2 =(d^n - x n) 2 

=(<9 M n) 2 - 2#<9 M n • x n) + x n) • x n) 

= (<V*) 2 - 2g^ ■ (n x d^n) + • ^)(n ■ n) - g 2 {^ ■ n) 2 

=g 2 Ml + 2g 2 (C^ + B M + X„) • B M + £ 2 (C M + B M + X M ) 2 - /(C^) 2 
=(? 2 lJ + 2^ 2 lJ + 2<?% • B M + g 2 (M ll + X M ) 2 

=VB 2 + 4</% ■ B M + g 2 Xl = g 2 {2M ll + X M ) 2 , (B.33) 

and there is no guarantee for the invariance of (2B M + X M ) 2 under the gauge transfor- 
mation II. 

Note that in the usual continuum limit the signature in front of ^ is not impor- 
tant. However, in this case, we can not adopt the 

U x ^ = exp{ieg^(x)} = exp{ze(7[C M (x) + B M (x) + X M (x)]}. (B.34) 

The identification (B.12) is necessary in order to recover the continuum expressions 
in the naive continuum limit as shown below. This is because the odd term in ^ 
plays the important role in this case and we can not change the signature arbitrarily. 

However, if we adopted (B.34), then the naive continuum limit left the B M depen- 
dent term after the CFN decomposition, 

=j ~ \e 2 {(d^n) 2 + 2g(d IM n x n) ■ ^ + g 2 {n ■ n){^ ■ sfj - g 2 {n • ^) 2 } + 0(e 3 ), 
= f -\e 2 m-^]n) 2 + 0(e 3 ) 

= y " ^V(2B M + X M ) 2 + 0(e 3 ). (B.35) 
Note that 

~ - \e 2 g 2 {K " O ~ (» • + 0(e% 

= f -^V(B M + X M ) 2 + 0( e 3 ). (B.36) 

If we wish to adopt (B.34) as the definition of the link variable, we must change the 
definition of B M as 

Bju = # _1 n x ^n, (B.37) 

which is determined from the condition of the covariant constant for the different 
co variant derivative, 

D^Vjn := d^n - g Y fl xn = d^n -jl^xn := [B] n = 0. (B.38) 
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Then we obtain the naive continuum limit, 

= y - ^e 2 {(^n) 2 + 2(7(9 M n x n) • ^ + # 2 (n • n)(^ ■ ^) - g 2 (n ■ ^f) + 0(e 3 ), 
=§ "^ + 0(4 (B.39) 
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